Large Mixing Angle Sterile Neutrinos and Pulsar Velocities 
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We investigate the momentum given to a protoneutron star, the pulsar kick, during the first 10 
seconds after temperature equilibrium is reached. Using a model with two sterile neutrinos obtained 
by fits to the MiniBoone and LSND experiments, which is consistent with a new global fit, there is 
a large mixing angle, and the effective volume for emission is calculated. Using formulations with 
neutrinos created by URCA processes in a strong magnetic field, so the lowest Landau level has a 
sizable probability, we find that with known parameters the asymmetric sterile neutrino emissivity 
might account for large pulsar kicks. 
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I. INTRODUCTION 

The existence of sterile neutrinos, in addition to 
three active neutrinos of the standard model, is of 
great interest for both particle physics and astro- 
physics. The Los Alamos LSND experiment found 
evidence for antincutrino oscillation lj. An analy- 
sis of LSND, along with other short-baseline experi- 
ments, foundQ that the standard model is not con- 
sistent with the data and the best model has two 
sterile neutrinos in addition to three active neutri- 
nos. 

The recent experiment by the MiniBooNE 
Collaboration^ found that the data for electron 
neutrino appearance showed an excess at low ener- 
gies, in comparison to what was expected in the stan- 
dard model. This data, along with the LSND data, 
has been analyzed in a model with two light sterile 
neutrinos [H, and compared to MiniBooNE data@. 
The mixing angles of two light sterile neutrinos were 
extracted. (See, however, Ref.@, which questions 
the accuracy of the r esults of RefH). An analy- 
sis of MiniBooNE, LSND, and other experiments 
based in part on earlier work with the seesaw 
mechanism [||], also found that the preferred model 
to fit the data consists of two light sterile neutrinos, 
in addition to three active neutrinos. 

Very recently a global analysis of updated Mini- 
BooNE data, along with data from LSND, KAR- 
MEN, NOMAN, Bugey, CHOOZE, CCFR84, and 
CDHS has been carried out@|. There are two pos- 
sible scenerios, both compatible with the results of 



Refs.Jl G|. Fits to appearance and disappearance 
with a 3+1 hypothesis have very good x 2 probabil- 
ities, and also find CPT violation. The CPT vio- 
lation, with different mixing angles and mass dif- 
ferences for neutrinos vs antineutrinos, might be 
understood as matter effects, similar to the MSW 
effect pil [U, but this has not yet been proved. Also 
a 3+2 model was shown to give a good \ 2 probability 
and is compatible within the range of the parameters 
for the 3+2 model used in Refs.0,JH], but without 
the uncertainty pointed out in Ref.[6(. 

In the present work we apply this range of fits[l,[|| 
to MiniBooNE /LSND to the study of pulsar veloc- 
ities. The gravitational collapse of a massive star 
often leads to the formation of a neutron star, a pul- 
sar. It has been observed that many pulsars move 
with linear velocities of 1000 km/s or greater, called 
pulsar kicks. See Ref.[l2| for a review. Due to the 
high density, the active neutrinos have a small mean 
free path, and only escape at the surface of the neu- 
trinosphere. During the first ten seconds after tem- 
perature equilibrium is attained during the gravita- 
tional collapse of a massive star, the main cooling is 
via emmision of neutrinos produced by the URCA 
process, and during the next ten seconds via neutri- 
nos produced by the modified URCA process. We 
have investigated the pulsar kicks which arise from 
the modified URCA processes in the time interval 
10-20 sec after the supernova collapses, when the 
neutrinosphere is just inside the protoneutron star. 
With a strong magnetic field and temperature, so 
that the population of the lowest Landau level is ap- 
proximately 0.4 of the total occupation probability, 
and we find large pulsar kicks [13j. 

The largest neutrino emission after the super- 
nova collapse takes place during the first 10 sec- 
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onds, with URCA processes dominant. The possi- 
bility of pulsar kicks from anisotropic neutrino emis- 
sion due to strong magnetic fields during this time 
was discussed more than two decades ago[T3]. It 
has been shown that, with the strength of the 
magnetic field expected during this period, the low- 
est Landau level has a sizable occupation probabil- 
ity, which produces the neutrino emission asymme- 
try that is needed for pulsar kicks. However, due 
to the high opacity for standard model neutrinos in 
the dense region within the neutrinosphere, few neu- 
trinos are emitted, and the large pulsar kick is not 
obtained [l6l|. 

Sterile neutrinos with a small mixing angle have 
small opacities. It has been shown [171] that, using 
the model of Ref. [l5[ and assuming the existence of 
a heavy sterile neutrino (mass > 1 kev) with a very 
small mixing angle constrained to fit dark matter, 
the pulsar kicks could be explained. More recently 
the effects of such sterile neutrinos with large masses 
and small mixing angles have been studied for other 
processes [l8j. See Ref.[l9T] for a review of processes 
that moght be associated with dark matter sterile 
neutrinos. 

In the present paper we use the fits of Refs.Jl, [|| 
with two sterile neutrinos to investigate the possi- 
bility of obtaining the large pulsar velocities which 
have been observed. The values of the masses and 
mixing angles are within the range found in calcula- 
tions based on the seesaw mechanismQ- Note that 
our model differs from that of Ref[17] in that with a 
much larger mixing angle there is a higher probabil- 
ity of sterile neutrinos, but a much smaller effective 
volume, due to a larger opacity. However, as we 
shall show, since the mean free path is much larger 
than those of standard neutrinos, under the condi- 
tions in which standard neutrinos produce a pulsar 
velocity of 2-300 km/s, the MiniBoone/LSND sterile 
neutrinos can give a kick of more than f 000 km/s. 



II. ASYMMETRIC STERILE NEUTRINO 
EMISSIVITY AND PULSAR KICKS IN 
LIGHT TWO-STERILE NEUTRINO MODEL 

Within about f second after the gravitational col- 
lapse of a large star, the neutrinosphere is formed 
with a radius of about 40 km, with temperature 
equilibrium. Within about 10 seconds about 98% of 
neutrino emission occurs, with neutrinos produced 
mainly by URCA processes. Due to the strong mag- 
netic field, neutrino momentum asymmetry is pro- 
duced within the neutrinosphere, but with a small 



mean free path they are emitted only from a small 
surface layer of the neutrinosphere, and the pulsar 
kick cannot be accounted for. If a standard active 
neutrino, say the electron neutrino, oscillates into a 
sterile neutrino, it will escape from the protoneu- 
trino star and neutrinosphere, unless it oscillates 
back into the active neutrino. The mixing angle 
plays a key role. In the work of Fuller et al[17| the 
mixing angle is so small that the sterile neutrinos 
are emitted. In the present work the starting point 
is the analysis of MiniBooNE and LSND data, with 
the two or more sterile neutrinos with small masses 
and large mixing angles. Before we can proceed, 
however, it is essential to determine possible effects 
of the high density and temperature of the medium 
on the mixing angles. 



A. Mixing Angle in Neutrinosphere Matter 

It has long been known that dense matter can af- 
fect neutrino states. The famous MSW effect fTol. fTTI] 
for understanding solar neutrinos, and the study of 
oscillations of high energy neutrinos [2(| are studies 
of mixing of active neutrinos in matter. There have 
been many other studies. In the present work we are 
dealing with sterile/active neutrino mixing given by 
the mixing angle 9 m in neutrinosphere matter 



\vi > = cos6 m \v e > —sinQ m \v s > 
\v 2 > = sin9 m \v e > +cos0 m \v s > 



(1) 



In the work Ref [171 it was shown that the mixing 
angle for sterile neutrinos that can account for dark 
matter as well as those produced in the neutron star 
core is almost the same as the vacuum value. Start- 
ing from the much larger mixing angles for the sterile 
neutrinos that seem to account for the MiniBoonE, 
LSND data, we need the value of the mixing angles 
in the neutrinosphere, as we discuss below. The ef- 
fective mixing angle in matter, 9 m can be related to 



the vacuum mixing angle, 9 by 21 



sin 2 (29 m ) 



sin 2 (29) 



sln ^29) + {cos{29)-^ 



(2) 



In Eq© V T is the finite temperature potential, 
while the finite density potential due to asym- 
metries in weakly interacting particles has been 
dropped as it vanishes when temperature equilib- 
rium is reached 21] . A convenient form for V T , with 
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the background of both neutrinos and electrons in- 
cluded, is given in Ref[23| 



V 1 = 



28ttGI 



F sin 2 9 w (l + 0.bcos z 8 w )pT i , (3) 



45a 



MiniBoone, LSND Data and Two Sterile Neutrinos 
M. Maltoni and T. Schwetz, hep-ph/0705.0107 (2007) 



with Gf,9w the standard weak interaction param- 
eters and a = 1/137. Assuming T=20a MeV, with 
a <1.0, p=6MeV, and {5m) 2 = 1.0 ev 2 0,[i|, we find 



Figure by M. Sterbcnz, LANL 



2pV 7 
(5m) 



5.1 x 1(T 3 6V < cos{29) . (4) 



Therefore, the mixing angle in the neutrinosphere 
medium is approximately the same as the vacuum 
mixing angle. This agrees with Ref(l7l|. 
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B. Emissivity With a Light Sterile Neutrino 



We now use the fits to MiniBooNe and LSND with 
light sterile neutrinos to estimate pulsar kicks. The 
MiniBooNE results are consistent with the LSND 
results and CPT only if there are at least two ster- 
ile neutrinos. Models with three sterile neutrinos 
have also been considered [y, la]- Fits to the Mini- 
BooNE experiment and the LSND results by Ref. [1] 
in Ref. [H with two sterile neutrinos are shown in Fig. 
1. 

From the Sterbenz/Maltoni-Schwetz fits one finds 
for the mixing angles of the two sterile neutrinos: 

(sin26 ls ) 2 = 0.004 

(sin29 2s ) 2 = 0.2, (5) 

and the masses are negligibly small. Note that this 
is in contrast to the parameters of Ref. [13], with 
the constraint of dark matter giving a mixing an- 
gle of (sin29d m ) 2 — 10~ 8 , and a mass greater than 
1 keV. In our present work we will use values for 
(sin29) 2 in the range 0.2 to 0.004 to estimate the 
pulsar kick, which is compatible with the recent 
global analysis @ 

The probability of asymmetric emission, giving a 
pulsar kick, does not depend directly on the sterile 
neutrino mass in our model, but is proportional to 
the (sin29 s ) 2 . It is the large mixing angles found in 
fits to MiniBoone and LSND that lead us to carry 
out the investigation in the present paper. 



FIG. 1: sm 2 (26»i s )=0.004; sm 2 (26» 2s )=0.2 

C. General Formulation of Neutrino 
Emissivity With a Strong Magetic Field 

The neutrino emissivity is give n in general form 
in many papers, e.g., see Refs [23l l24| : 

„ n 4 f d 3 p* d 3 q» f d 3 q e 
1=1 J (27r) 3 2^(27r) 3 J (2tt) 3 

( 27r ) 4 E ^™ iM ( 6 ) 

5(Efi na [ Ei n iti a i)5(pfi na i Pinitial) ? 

where M is the matrix element for the URCA pro- 
cess and T is the product of the initial and final 
Fermi-Dirac functions corresponding to the temper- 
ature and density of the medium. The main source 
of the asymetric emissivity that produces the pul- 
sar velocity is the fact that the electron has a large 
probability to be in the lowest (n=0) Landau level. 
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See Refs [25|, |26[ for a discussion of Landau levels. 
The asymetric emissivity can be seen by considering 
the weak axial intereaction, Wa, 



G - 
I„ = *( g e ) 7M (l- 75 )*(^) 



(7) 



with G = 3a = 1.26, the x are the nu- 

cleon spinors, and the lepton wave functions are 
^>(q e ), 1 i'(q u ), where q e and q v are the electron and 
antincutrino momenta, respectively. The key to the 
asymmetric emission is given be the trace over the 
leptonic currents, Tr[ljlj], 



d 2 qlTr[ltlj] ~ 8irE e [( q y 5 i3 + (q»y6 j3 



-W?W = B = z), (8) 

with the magnetic field B in the z direction. We only 
consider the weak axial force, which is dominant. 
Using the relationship given in Eq([8|) one can show 
that the result of the traces and integrals over the 
axial product matrix element has the form (P = z) 



// |M ' 



oc 



(9) 



Details are given in Ref [l3|, where it is shown that 
the asymmetric neutrino emissivity, using the gen- 
eral formulation of Ref 1231 , is 



0.64 x 10 21 T 9 7 P(0) x / 



erg cm 3 s 1 



1 



Pn 



V ef fAt 



(10) 



where Tg = T/(1Q 9 K), p ns is the neutron star mo- 
mentum, P(0) is the probability of the electron pro- 
duced with the antineutrino being in the lowest Lan- 
dau state, f=.52 is the probability of the neutrino 
being at the + z neutrinosphere surface [13] ■ V e ff 
is the volume at the surface of the neutrinosphere 
from which neutrinos are emitted, and At ~ 10s is 
the time interval for the emission. 

We derive P(0) and V e ff, the effective volume for 
the emissivity, in the next two subsections. 



D. P(0) = Probability for the Electron to be in 
the n=0 Landau Level 



Just as in our previous work in which Landau 
levels play a crucial roleJTJj], only the lowest Lan- 
dau level, for which the helicity is -1/2 (rather than 



±1/2 as with the usual Dirac spinors) gives asym- 
metric emission. The probability that an electron in 
a strong magnetic field is in the lowest (n=0) Landau 
level, P(0), can be calculated from the temperature, 
T, and the energy spectrum of Landau levels [2o1.[26T|. 
A particle with momentum p and effective mass m* 
in a magnetic field B in the nth Landau level has the 
energy 



E L {p,n) 



(m* 



) 2 + 2(m e *) 2 ^-n, (11) 



with B c = 4 x 10 13 G, and m* is the effective mass of 
the electron at the high density of the protoneutron 
star and neutrinosphere. 

From standard thermodynamics the probability of 
occupation of the n=0 Landau level, P(0), is given 
by [H: 



P(0) 



F(0) 



F(0)+2£~F(n) 



(12) 



where F(n), with magnetic field B, temperature T, 
and chemical potential fi, is 



F(n) 



dp 



[m n -m p - E L (p, n)} 2 
l + exp[(E L (p, n)~fi)/T] 



(13) 



The electron energy is restricted to magnitudes 
greater than /j,, but the integrals in Eq.([7]) are in- 
sensitive to Pmin, so we take p m i n =0 as in Ref[l7l|. 

We agree with the estimate of Ref. [17] for P(0). 
Note that if we had used the free electron mass, m e , 
in the Landau energies (Eq.([5])) we would have ob- 
tained a much smaller value for P(0). For B= 10 16 
G, n = 40 MeV, m*=4 MeV and T u - Bphere = 20 
MeV, P(0) ~ 0.3. This is similar to our estimate 
of P(n=0) ~ 0.4 at the surface of the protoneu- 
tron star at about 10 seconds [l3|. Therefore our 
result for asymmetric emissivity differs from that of 
Fuller et al [17| mainly in that we have a much larger 
mixing angle, and a much smaller effective volume, 
since the sterile neutrinos oscillate back to active 
neutrinos within the neutrinosphere; and therefore 
our emission only takes place near the surface of the 
neutrinoshere. However, in contrast with purely ac- 
tive neutrino emission in which the opacity results 
in very small pulsar kicks [l6j |. the sterile neutrinos 
have a much larger effective volume, and can there- 
fore produce much larger pulsar velocities 
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E. Estimate of V e ff— Effective Volume for 
Emission 

To estimate V e ff we make use of the early study 
of opacity in about the first 20s of the creation of 
a neutron star via a supernova collapse [U [H, [3(| , 
and a recent detailed study of neutrino mean free 
paths [28j| . Since the mean free path of the sterile 
neutrino is determined by that of the standard neu- 
trino to which it oscillates, A, we make use of studies 
of active neutrino mean free paths. First note that 
the neutrino mean free path is given by 



1/A = 



d 3 p 



M fi [l-n(q)] 



(2tt) 3 

x{l + e ^-E v )/kT ) 



(14) 



where Mfi is the weak matrix element and n(q) is 
the Fermi distribution. For the calculation of the 
sterile neutrino, for which Mfi = 0, one can use 
the value of 1/A with a factor of sin 2 (20) from the 
matrix element and another such factor from the 
occupation probability. From the results of previous 
authors, for T in the 10 to 20 Mev range and u in 
the 20 to 40MeV range, we estimate that A ~ 1.0cm 
This gives a range for the effective sterile neutrino 
mean free path 



A., 



5.0 to 250 cm 



(15) 



For a neutrinosphere radius of 40 km, with A s < < 
R v this gives us V eff = (4tt/3)(J£ - (R v - A,) 3 ) ~ 

From Eq.fTtj]). R v = 40 km, and A=1.0 cm, 
067x10^ cm 

p ns = M ns V ns ~ ■ 2(nf)\ T 9 9 m j (16) 

with Tg = lt J) K ■ Taking the mass of the neutron 
star to equal the mass of our sun, M ns = 2 x 10 33 
gm, we obtain for the velocity of the neutron star 



For example, for T=10 Mev =1.16 X 10 11 kK, 

km km 

v ns ~ 47.3 to 2, 370 , (18) 

s s 

which means that sterile neutrino emission could ac- 
count for the large pulsar kick with the parameters 
extracted from Refs 0, [|[. If we use the physical 
parameters that give Eq. p^)) for electron neutrinos, 



we obtain a pulsar velocity of v ns = 95 km/s, which 
is consistent with previous predictions by several au- 
thors. 

It should be noted that the study of the Mini- 
BooNE and LSND results are in progress, and the 
mixing angles that result could be different from 
those obtained in Refs. [I, For this reason wc 
have used a range of parameters. Preliminary data 
from the MiniBooNE/Minos experiment [3l|, how- 
ever, is consistent with the MiniBooNE results [31. 
We also once more point out that although Ref Q 
questions the accuracy of the peramaters extracted 
by Ref our model is compatible with the recent 
global analysisQ, with very good \ 2 probabilities. 



III. CONCLUSIONS 

Because of the strong magnetic fields in protoneu- 
tron stars and the associated neutrinosphere, the 
electrons produced in the URCA processes that 
dominate neutrino production in the first 10 seconds 
have a sizable probability, P(0), to be in the low- 
est (n=0) Landau level. This leads to asymmetric 
neutrino momentum. With the mixing angles found 
m Refs 4, 5], we find that the sterile neutrinos pro- 
duced during this period for high luminosity pulsars 
can give the pulsars velocities of greater than 1000 
km/s, as observed. We emphasize that our results 
are consistent sterile neutrino parameters based on 
a global analysis of data from MiniBoonE, LSND, 
KARMEN, NOMAN, Bugey, CHOOZE, CCFR84, 
and CDHS@|, rather than a model. There is a high 
probability that light sterile neutrinos with a large 
mixing angle exist, which is the basis of our work. 

There is a strong correlation of the pulsar velocity 
with temperature, T. Since it is difficult to determine 
T accurately, it is difficult for us to predict the ve- 
locity of a pulsar whose kick arises from sterile neu- 
trino emission. On the other hand, if the pulsar kick 
arises from the asymmetric emission of active neutri- 
nos produced by the modified URCA processes after 
10 seconds, also proportional to P(0)[1J|, then T can 
be determined by an accurate measurement of the 
neutrinos from the supernova. Therefore, in future 
years, with much more accurate neutrino detectors, 
one could predict the velocity of the resulting pulsar. 
Unfortunately, the energy of emitted sterile neutri- 
nos cannot be measured. From our results in the 
present paper and those in Ref[13[, high luminosity 
pulsars receive a large kick both from sterile neutri- 
nos in the first ten seconds and standard neutrinos 
in the second ten seconds. 
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